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ABSTRACT 

Natural  hazards  occupy  the  essential  and  regional  levels,  hence,  they  are  raised  as  a priority  issues.  The  2009 
Saudi  Arabia  floods  affected  Jeddah,  on  the  red  sea  (western)  coast.  As  of  January  3rd,  2010,  122  people  are  reported  to 
have  been  killed.  Roads  were  under  a meter  of  water.  Unfortunately,  Lack  of  knowledge  in  water  flow  modeling 
contributes  to  prevent  manage  flood  risks.  This  paper  provided  the  spreadsheet  model  as  solver  to  solve  wave  propagation 
due  to  sudden  closure  of  the  downstream  gate  or  closed  road.  A Wave  propagation  flow  is  simulated  by  numerically 
solving  the  one-dimensional  saint  venant  equations  by  using  a second-order  explicit  finite -difference  (McCormack) 
scheme.  This  method  has  been  verified  with  using  lux  diffusive  scheme.  Spreadsheets  as  new  solution  technique  is  used  to 
simulate  unsteady  flow  in  open  channel  which  is  a practical  method.  The  results  obtained  by  using  the  Saint  Venant 
equations  are  compared  to  determine  the  depth  at  each  elven  sections  which  the  maximum  depth  was  6.858437  m at  time 
around  1050  second.  It  is  found  that,  the  Saint  Venant  equations  give  sufficiently  accurate  results  for  the  maximum  flow 
depth  and  the  time  to  reach  this  value  at  a location  downstream. 

KEYWORDS:  Saint  Venant,  Spreadsheets,  McCormack  Scheme 

INTRODUCTION 

Transient  or  unsteady  flow  in  open  channels  can  be  called  propagation  of  water  waves.  Mathematically,  water 
waves  propagation  can  be  described  by  the  Saint- Venant  equations  (assuming  hydrostatic  pressure  distribution  in  the 
vertical  direction)  or  by  other  equations  such  as  simplified  wave  approximations,  non-inertia  wave,  gravity  wave  and 
dynamic  wave  models.  Since  the  1970s  different  methods  of  numerical  approaches  have  been  developed  in  order  to  solve 
Saint- Venant  equations  or  others,  such  as  Liggett  and  Woolhiser  (167),  Wylie  (1970)  and  Woolhiser  (1977).  The  analysis 
of  propagation  of  water  waves  is  vital  for  maintenance,  operation  and  emergency  planning.  Water  wave  propagation  is 
studied  by  numerical  solution  of  the  Saint- Venent  equations.  This  paper  will  focus  on  numerical  solutions  used  for  solving 
the  Saint-Venent  equations  that  describe  and  simulate  water  waves,  water  wave  propagation  and  dam  break  simulation. 
The  main  goal  of  this  paper  is  to  understand  water  wave  propagation. 

Computational  hydraulics  is  an  important  for  water  resource  engineers,  it  is  a tool  of  simulation  water  flow  under 
certain  condition  to  ease  studying  and  visualize  a certain  case  . In  this  paper  McCormack  scheme  is  used  for  simulating 
sudden  closure  of  downstream  gate  and  lux  diffusive  scheme  in  order  to  verification  .recently  some  of  published  paper 
used  spreadsheet  (Excel)  in  solving  complex  3 dimension  numerical  problem  particularly  in  hydraulics  field  for  instance 
HALIL  KARAHAN,  M.  TAMER  AYVAZ  utilize  spread  cheat  (excel)  in  Groundwater  Modeling  .spread  cheat  (excel) 
flexible  and  workable  to  use  which  is  showing  the  result  numerically  and  graphically.  All  numerical  methods  are  subjected 
to  stability  constraints  which  restrict  the  values  allowed  for  the  time  step  for  a given  grid  [1,2].  For  explicit  schemes,  this 
can  result  computationally  long  .Better  explicit  schemes  or  implicit  schemes  overcome  this  disadvantage  but  with  more 
complexity  [2]. 
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The  saint  venant  equations  or  dynamic  equation  appear  in  many  forms  in  the  literature  and  can  either  be  written  as 
a set  of  integral  or  differential  equations  [4],  the  Saint  Venant  equations  are  a special  case  of  the  Navier-Stokes  equations. 
In  case  of  one  dimension  unsteady  flow,  momentum  equation  or  energy  equation  and  continuity  equation(PDE)  must  be 
solved  which  are  having  at  least  2 unknown  variable  the  following  method  are  used  for  solving  them  . 

• Characteristic  method 

• Kinematic  wave  simulation 

• Diffusive  wave  simulation 

• Full  dynamic  simulation 

• Finite  deference 

• Finite  element 

• Finite  volume 

• Finite  boundary 

PROBLEM  STATEMENT 

The  propagation  of  water  waves  might  cause  dam  failure  and  can  cause  significant  damage,  particularly  when 
dams  are  located  in  a residential  area,  such  as  is  the  case  in  Jeddah.  Water  wave  propagation  is  not  completely  understood, 
and  from  this  point  of  view  it  is  important  to  induce  researchers  to  study  this  phenomenon. 

METHODOLOGY 

Seeking  the  scientific  approach  of  using  numerical  solution.  Deductive  and  inductive  is  adapted  in  this  research 

OBJECTIVE 

This  paper  aims  to  improve  the  level  of  understanding  of  the  propagation  of  water  waves  in  open  channels.  The 
main  approach  of  this  paper  will  be  numerical  solutions  used  for  solving  the  Saint -Venant  equations  that  describe  and 
simulate  water  waves,  particularly,  transient  condition  which  are  produced  by  the  sudden  closure  of  a downstream  gate  or 
closed  road.  Using  spread  cheat  (excel)  to  simulate  water  waves. 

FUNDAMENTAL  EQUATION  OF  UNSTEADY  FLOW 

The  flow  of  the  water  can  be  described  mathematically  based  on  three  fundamental  laws  of  physics  as  follows: 

• Conservation  of  matter  (or  mass). 

• Conservation  of  energy,  (Bernoulli  equation) 

• Conservation  of  momentum. 

The  concepts  is  described  briefly  and  given  below. 

GOVERNING  EQUATIONS 

Continuity  equation  and  momentum  or  energy  equation  are  the  common  governing  equations  which  together  form 
the  differential  version  of  the  saint  venent  equations  see  equations.  As  our  case  assumption  that  the  system  has  no  energy 
losses,  the  covering  equations  are  Continuity  equation  and  momentum  equations.  By  solving  them  we  will  obtain  the 
variable  in  term  of  each  other. 


Numerical  Solution  of  the  Shallow  Water  Equation  by  Using  Spreadsheet 
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Momentum  or  Dynamic  Equation 

By  applying  Newton’s  2nd  law  to  our  elemental  length  of  channel  we  have  f=mg.  In  case  of  unsteady  flow  here  is 
the  equation: 


6Q  dQV  Acy 

— - + — — + gA— 
dt  8x  8x 


gA(s0-sf) 


(1) 


Substitution  of  Q = AV  into  Eq.  1,  equation  leads  to: 


8 V 
dt 


+ V 


8V 

6x 


+ = s(s,-s,) 


Figure  1:  Definition  Sketch  for  Derivation  of  St.  Venant  Equations 


Continuity  Equation 

Unsteady  flow  can  be  written  depends  on  mass  (matter)  conservation  as  flows: 

Substitution  of  Q = AV  into  equation,  Continuity  Equation  leads  to  be: 

3A  5AV  „ 

1 = o 

ar 

One  dimensional  continuity  equation  can  be  written  as  (Chaudhry  2002): 

3 t \ /T  J o'jc  a'x 


(2) 


Equation  for  celerity  of  a small  gravity  wave  in  a trapezoidal  channel  can  be  written  as: 

C = V g*T/A 

A FIRST  NUMERICAL  METHOD  OF  SOLUTION 

Eqs.  (1)  and  (2)  are  nonlinear  analytical  solutions  which  are  both  partial  differential  equations,  a finite-difference 
method  is  used  to  solve  these  two  governing  equations.  Therefore,  a numerical  approach  is  used  here  for  their  solution.  It  is 
necessary  to  use  third  or  higher-order  accurate  numerical  methods  to  solve  these  equations  (Abbott  1979).  In  this  paper,  the 
method  developed  by  McCormack  and  lux  (1976)  (Chaudhry  2003)  are  used  for  solving  the  saint  venet  equations  for 
studying  the  wave  of  downstream  dam  gate  sudden  closure. 
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Two  well-known  numerical  methods  will  be  used  in  this  case  although  they  are  both  approximations.  One  is 
McCormack  sachem  and  the  anther  is  lux  diffusive  sachem  which  is  both  based  on  using  the  ideas  of  where  information 
coming  from  which  itself  determined  by  the  characteristics  see  figure  [2]. 

Consider  the  arrangement  of  grid  shown  in  Figure  3,  in  which  the  nodes  are  spaced  at  regular  Ax  intervals  along 
the  direction  of  the  x axis  and  at  At  intervals  along  the  t axis  [2]. 


First,  an  intermediate  flow  field  (v  velocity  ,h  water  depth)  is  computed  by  a predictor-corrector  procedure  see 
appendix  , The  partial  derivatives  of  the  dependent  variables  can  be  written  as: 


Left  boundary  conditions 


*’*-0  * yrb  - 
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Right  boundary  conditions 
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The  left  and  right  boundry  condtions  are  shown  on  the  figure  2. interior  value  of  v,  y for  the  new  time  step  can  be 
written  as  predictor  and  then  can  be  corrected  to  be  called  predictor  , after  that  , by  taking  the  average  between  the 
corrector  and  predictor  we  get  the  target  value  , the  posture  as  flows: 

Predictor 

Substitution  of  these  expressions  for  the  partial  derivatives  and  simplification  and  predictor  can  be  written  as 


follows: 


Numerical  Solution  of  the  Shallow  Water  Equation  by  Using  Spreadsheet 
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u*  = JJf  - 


At 

Ax 


(F?  - FjLx)  - S’?  At 


vr  = V;"-  v>-£* y^t  {So- % 


(3) 


Corrector 

Substitution  of  these  expressions  for  the  partial  derivatives  and  simplification  and  Corrector  can  be  written  as 


follows: 


, /At 

u?*  = u*  - 2^ 


(FJ+1  - F *)  - SJZit 


vr,  (4-  sfi ) 


Where  , 


„ a*  . f , d’ 


T 


Finally  , interior  vlue  of  v,y  can  be  estimate  by  taking  the  average  of  predictor  and  corector  as  fllows  : 


U 


H-  XJ**) 


A SECOND  NUMERICAL  METHOD  OF  SOLUTION 


(4) 


(5) 


As  mention  before  lux  schemes  is  used  to  verify  McCormack  sachem,  it  is  a method  of  approximation,  lux 
schemes  has  same  boundary  conditions  which  are  the  equation  of  positive  and  negative  charstarstic  line  see  figure  2, 
whereas  interior  point  of  y,  v can  be  written  as  ( chaudhry),  2003): 


y^x 


= +2/ti)  - 

~ 2 A'X  (y*-+- * ~ yi~  1 ) 

v**~*~1  = 2 + ^+1)  ~ k ^~y(y^x  — y?- 1) 


3.  At 
2 


2 Ax' 


K*(KXi  — K-i)  -h  gAt(S0  — S}) 


(6) 


Where, 
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STABILITY  AND  ACCURACY 

Stability  of  both  schemes  is  controlled  by  At  and  Ax  as  follows: 

• The  positive  characteristic  through  L has  slope  dt/dx  = l/(v+c). 

• The  negative  characteristic  through  R has  slope  dt/dx  = l/(v-c). 

So,  to  keep  within  this  domain,  these  criteria  for  the  positive  characteristic  must  be  meeting:  dt/dx  > (v+c) 


Figure  3:  Characteristics  around  Point  p 
SIMULATION  BY  USING  SPREADSHEETS 

Solving  the  PDEs  by  using  spreadsheet  is  one  of  the  best  technic.  Graphical  interface,  numerical  and  visual  result 
and  fast  calculating  capabilities  are  one  of  the  spreadsheets  important  advantages. 

Changing  any  input  will  be  festally  reflected  to  the  graph.  FDM  equation  can  be  written  easily  then  copied  to 
other  cells  without  writing  the  equations  again  to  all  cells  separately  which  represent  the  iteration  in  spreadsheets.  The 
objective  of  using  spreadsheets  in  this  paper  is  to  use  flexible  tool  in  order  to  water  flow  Simulation. 

The  following  input  parameters  are  used  for  trapezoidal  are:  channel  length  from  upstream  to  the  dam, 
L = 5000  m;  dam  location,  L = 5000m;  initial  water  depth  in  channel  y = 5.79m;  acceleration  due  to  gravity  g=9.81  m2/s; 
Manning  roughness  coefficient,  n=  0.013,  bed  slope,  So  = 0.00008;  and  step  size,  dx  = 500m;  channel  width,  6.1m; 
channel  lateral  slope,  z=1.5;  initial  steady  state  discharge,  Q=126. 

In  order  to  obtain  a stable  solution  for  explicit  simulation  can  be  seen  in  cell  F 4,  5 in  figure  4,5  where  dt  is  equal 
to  67  according  to  the  stability  condition  shown  above  . 

The  solution  domain  for  the  wave  propagation  due  to  sudden  closure  of  the  downstream  gate  based  on  lux 
scheme,  which  uses  FDM,  can  be  seen  in  Figure  4. 


Numerical  Solution  of  the  Shallow  Water  Equation  by  Using  Spreadsheet 
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Figure  4:  Solution  Domain  of  the  Lux  Scheme 

Depending  on  our  requirements  we  can  easily  generate  FDM  equations  in  the  cell  as  much  as  need  based  on  the 
Dx  and  Dt,  which  are  the  size  of  grids  in  the  solution  domain  see  attached  spreadsheet. 

For  lux  scheme,  Solution  of  FDM  equations  is  carried  out  based  on  Eq.  (6)  in  the  following  spreadsheet  format: 

Left  Boundary  Condition  Excel  Code 

(velocity)E27=F20-((SQRT($C$12*F17/F18))*F15)+($C$12*$F$4*($C$6 

F21))+((SQRT($C$12*F17/F18))*E22) 

Y=  constant  (5.79) 

Right  Boundary  Condition  Excel  Code 


V=0 


(flow  depth)  022  =(((SQRT($C$12*N17/N18))*N15)+($C$12*F4*(C6-N21))+N20)/(SQRT($C$12*N17/N18)) 

Interior  Value  of  V,  Y Excel  Code 

(velocity)I27=0.5*(H20+J20)-(0.5*$F$5*$C$12*(J15-H15))-(0.5*$F$5*0.5*(J20+H20)*(J20- 

H20))+($C$12*$F$4*($C$6-(0.5*(J21+H21)))) 

(flow  depth)  =0.5*(H15+J15)-(0.5*$F$5*0.5*(H19+J19)*(J20-H20))-(0.5*$F$5*0.5*(J20+H20)*(J15-H15)). 

For  McCormack  sachem,  Solution  of  FDM  equations  is  carried  out  based  on  Eq.  (3,4,5)  in  the  following 
spreadsheet  format: 

Left  Boundary  Condition  Excel  Code 

(predictor  velocity)  E34=F24-((SQRT($C$12*F19/F20))*F15)+($C$12*$F$4*($C$6- 
F26))+((SQRT($C$  1 2*F1 9/F20))*E27). 
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(corrector  velocity)  E35=E24-($F$5*(E34)*(F34-E34))-($C$12*$F$5*(F28-E28))-($C$12*$F$4*($C$6-E25)) 
Y=  constant  (5.79). 

Right  Boundary  Conditions  Excel  Code 


V=0 


(predictor  flow  depth)  028  =015-($F$5*021*(024-N24))-($F$5*024*(015-N15)) 

(corrector  flow  depth)  029  =(((SQRT($C$12*N19/N20))*N15)+($C$12*$F$4*($C$6- 
N25))+N24)/(SQRT($C$12*N19/N20)). 

Interior  Value  of  V,  Y Excel  Code 

(predictor  velocity)  H34  =H24-($F$5*(H24)*(H24-G24))-($C$12*$F$5*(H15-G15))-($C$12*$F$4*($C$6- 
H26)). 

(corrector  velocity  ) H35  =H24-($F$5*(H34)*(I34-H34))-($C$12*$F$5*(I28-H28))-($C$12*$F$4*($C$6-H25)) 
(predictor  flow  depth)  H28  =H15-($F$5*H21*(H24-G24)M$F$5*H24*(H15-G15)). 

(corrector  flow  depth)  H29  =H15-($F$5*H30*(I34-H34))-($F$5*H34*(I28-H28)). 

Average  velocity  H26  =(H34+H35)/2 
Average  flow  depth  H27=(H28+H29)/2 
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Figure  5:  Solution  Domain  of  the  McCormack  Scheme 

One  -dimensional  representation  of  hydraulic  heads  can  be  seen  in  Figure  6 for  McCormack  sachem  and  Figure  7 
for  lux  scheme.  This  figure  shows  the  level  of  hydraulic  heads  over  the  solution  domain.  Each  point  of  Figure  6,  7 indicates 
the  different  levels  of  hydraulic  heads.  The  flow  head  of  the  wave  propagation  can  be  seen  in  figure  8 for  1050  second  at 
which  the  water  rises  to  reach  the  maximum  flow  head.  Figure  8 shows  the  variation  between  flow  depths  for  1050  second 
by  using  McCormack  sachem  and  lux  scheme. 
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Figure  6:  One-Dimensional  Variation  of  Flow  Depth  Level  (t=1050  Second)  by  Using  McCormack  Scheme 
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Figure  7:  One-Dimensional  Variation  of  Flow  Depth  Level  (t=1050  Second)  by  Using  Lux  Scheme 


Figure  8:  Variation  between  McCormack  and  Lux  Scheme 
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CONCLUSIONS 

Natural  hazards  occupy  the  essential  and  regional  levels,  hence,  they  are  raised  as  a priority  issues  in  road  design 
and  protection.  The  2009  saudi  Arabia  floods  affected  Jeddah,  on  the  red  sea  (western)  coast.  As  of  January  3rd,  2010,  122 
people  are  reported  to  have  been  killed.  Roads  were  under  a meter  of  water.  Unfortunately,  Lack  of  knowledge  in  water 
flow  modeling  contributes  to  prevent  manage  flood  risks.  This  paper  provided  the  spreadsheet  as  solver  used  for  solving 
finite  deference.  Finite  deference  of  boundary  conditions  and  interior  value  were  written  by  using  two  different  schemes 
(McCormack  and  lux  scheme).  The  calculations  have  been  done  iteratively.  One  of  the  main  purposes  of  this  paper  is  to  try 
to  match  or  meet  the  scientific  methodology  for  studying  such  phenomena  through  using  one  of  numerical  solution  scheme 
and  anther  for  verification.  The  variation  between  tow  schemes  is  quite  large.  Lux  scheme  has  been  compared  with  Hanif 
Chaudhry  results  and  found  that  the  result  is  same.  Figures  show  lux  scheme  is  more  logical  than  McCormack  scheme. 
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